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This article presents a modified free-piston Stirling heat engine configuration in which a thermal buffer 
tube is added to sandwich between the hot and cold heat exchangers. Such a modified configuration may 
lead to an easier fabrication and lighter weight of a free piston. To analyze the thermodynamic perfor¬ 
mance of the modified free piston Stirling heat engine, thermoacoustic theory is used. In the thermoacou¬ 
stic modelling, the regenerator, the free piston, and the thermal buffer tube are given at first. Then, based 
on linear thermoacoustic network theory, the thermal and thermodynamic networks are presented to 
characterize acoustic pressure and volume flow rate distributions at different interfaces, and the global 
performance such as the power output, the heat input and the thermal efficiency. A free piston Stirling 
heat engine with several hundreds of watts mechanical power output is selected as an example. The typ¬ 
ical operating and structure parameters are as follows: frequency around 50 Hz, mean pressure around 
3.0 MPa, and a diameter of free piston around 50 mm. From the analysis, it was found that the modified 
free-piston Stirling heat engine has almost the same thermodynamic performance as the original design, 
which indicates that the modified configuration is worthy to develop in future because of its mechanical 
simplicity and reliability. 

© 2011 Elsevier Ltd. All rights reserved. 


1. Introduction 

The Stirling engine is driven by the cyclic hot expansion and 
cold contraction of the working gas. It is one of the most efficient 
heat engines. The Stirling engine is an external combustion engine, 
and many kinds of heat sources could be used such as solar energy, 
industrial waste heat and biomass energy. In addition, the Stirling 
engines have advantages of smooth operation, low noise, and envi¬ 
ronmental friendliness [1]. 

A free piston Stirling engine (FPSE) is a variable form of Stirling 
machines by eliminating the crank-shaft device of the traditional 
Stirling engine. Its thermodynamic performance is greatly influ¬ 
enced by the moving mass and spring stiffness of the displacer 
assembly (i.e. free piston) [2]. Compared with the traditional Stir¬ 
ling engine, it could be self-started and has better reliability, higher 
efficiency as well as lower vibration. Therefore, the free piston Stir¬ 
ling engine is getting more and more attentions and might have 
many applications [3-6]. 

However, there are some drawbacks for the current FPSE. A 
large part of the displacer works in the high temperature region, 
which leads to several disadvantages including difficult manufac¬ 
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turing, assembling and sealing, etc. In order to overcome these dis¬ 
advantages, a modified FPSE with a thermal buffer tube is 
proposed. In this new configuration, the high-temperature part of 
the displacer is replaced by a thermal buffer tube. Meanwhile, sev¬ 
eral layers of anti-radiation stainless steel meshes are stacked in 
the thermal buffer tube to reduce the radiation heat loss. In addi¬ 
tion, an ambient heat exchanger is set at the ambient region of 
the displacer to make sure that the displacer operate at ambient 
temperature. This configuration can effectively reduce the diffi¬ 
culty of manufacturing, assembling and sealing of the displacer. 
Also, the lifetime of the displacer would be greatly improved. 

2. Modified free piston Stirling engine 

Fig. 1 shows the schematic of a conventional FPSE. The main 
components of the FPSE include: a hot end heat exchanger, a 
regenerator, an ambient heat exchanger, a displacer, a power pis¬ 
ton, and a linear alternator. As the engine works, part of the acous¬ 
tic power from the expansion space is used to drive the mechanical 
device. The other acoustic power enters the ambient end of the 
regenerator, and then amplifies. 

Fig. 2 shows the schematic of a modified FPSE. It uses a thermal 
buffer tube to replace the high-temperature part of the displacer. 
Its working principle is the same as a conventional engine. The 
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Nomenclature 



English letters 

Greek letters 

i 

v^l 

T 

temperature ratio 

X 

displacement amplitude (m) 

Po 

dynamic viscosity (kg/m s) 

M 

moving mass (kg) 

6 

phase (°) 

r h 

hydraulic radius (m) 

y 

adiabatic index 

S 

cross area of regenerator (m 2 ) 

p 

thermo efficiency 

T h 

heating temperature (I<) 

4> 

porosity 

V 

volume of thermo buffer tube (m 3 ) 

CO 

angular frequency 

d 

diameter of thermo buffer tube (m) 



k 

spring stiffness (N/m) 

Superscript 

l 

length of regenerator (m) 

* 

conjugate 

P 

amplitude of pressure wave (MPa) 



Pm 

damping coefficient (kg/s) 

Subscript 


To 

ambient temperature (I<) 

i(l-4) 

status points 

u 

amplitude of volume flow rate wave (m 3 /s) 



w 

acoustic power (W) 





Fig. 1. Schematic of conventional FPSE (1. hot end heat exchanger, 2. regenerator, 3. ambient heat exchanger, 4. linear alternator, 5. planar spring, 6. back space, 7. expansion 
space, 8. displacer, 9. compression space, 10. power piston). 



Fig. 2. Schematic of modified FPSE (1. hot heat exchanger, 2. regenerator, 3. ambient heat exchanger, 4. linear alternator, 5. planar spring, 6. back space, 7. expansion space, 8. 
thermal buffer tube, 9. mesh, 10. ambient releaser, 11. displacer, 12. compression space, 13. power piston). 


performance of the displacer is determined by the diameter of the 
rod, the spring stiffness of the planar flexure and the damping coef¬ 
ficient. The thermal buffer tube is a circular tube. The working gas 
in the thermal buffer tube acts as a long gas piston, transporting 
the acoustic power in forms of pressure and velocity waves from 
the hot region to the ambient region. 


However, there are several problems emerging with the intro¬ 
duction of the thermal buffer tube. First, the radiation heat loss 
from the hot region to the ambient region will increase signifi¬ 
cantly. Worse more, the gas convection occurring in the thermal 
buffer tube would deteriorate the heat loss therein. To reduce 
the heat loss, several layers of stainless steel meshes are placed 
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in the thermal buffer to decrease the radiation heat loss and make 
the flow uniform. Meanwhile, an ambient releaser is introduced to 
make sure that the displacer work in the ambient temperature. 


3. Mathematical model 


For all kind of oscillating flow heat engines, the dynamic feature 
plays an important role. In the following analysis, assuming that 
the working gas is helium and the heat transfer process is ideal. 
Then, the dynamic models of the regenerator, the thermal buffer 
tube, and the displacer are established, respectively. Since that 
attention is focused on how the heat is converted into acoustic 
power efficiently in the thermodynamic loop of the engine, the 
coupling between the engine and the linear alternator is not con¬ 
sidered in this work. 

Fig. 3 shows a simplified flow chart of the modified FPSE with a 
thermal buffer tube. The pressure wave and the volume flow rate 
wave at the inlet of the regenerator are (pi, U i), the parameters 
of the points 2-5 have the same symbols as point 1 except for 
the subscripts. 

It is known that from the momentum and mass balances: 


I Pi = P4 = P^ 

\ U x = U 4 - U 5 

The dynamic model of the regenerator [7] is: 
fp 2 \ _ (FPP\ ppm \ /pi \ 

V& 2 y VWP1 Fum){u,j 

where 

FPP1 = 1 + b). FPU 1 = - R °( T + 1 ) /( T , b) 
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The dynamic model of the thermal buffer tube [8] is: 


/ p3 \ _ ( FPP2 FPU2 \(P2\ 
\U 3 )-\FUP2 FUU2j[u 2 ) 

where 

FPP2 = 1, FPU2 = 0 


FUP2 = -i —, FUU2 = 1 

yp 0 

The dynamic model of the displacer is: 


/ p 4 \ / FPP3 FPU3 \ f P 3 \ 
\U 4 J ~ \FUP3 FUU3j\U 3 ) 

where 


FPP3 


A 3 

A~4 


FPU 3 = 


1 

2 ^ 32^4 


R m + ico 






FUP3 = 0, FUU3 = !K A 

^3 


A 3 and A 4 represent the area of the inlet and outlet of the displacer, 
respectively. 

From Eqs. (2) and (3): 


/ p 4 \ / FPP FPU\fpA 

[u 4 )-\fup fuu){u,J 



From Eqs. (l)-(5): 

( n — 1 -FPP f\ 

— FPU Pi 

U 2 = FUP\ * P! + FUU\ * U 1 
p 2 = FPP\ * p\ +FPm * U 1 
U 3 = FUP2 * p 2 + FUU2 * U 2 
p 3 = FPP2 * p2 + FPU2 * U 2 
U 4 = FUP*p 1 +FUU*U 1 
U 5 = U 4 - U\ 

, p 5 = p 4 = Pi 


From Eq. (6), it indicates that once the pressure wave of the in¬ 
let of the regenerator is known, the remaining parameters could be 
calculated easily. Then, the acoustic power and the thermal effi¬ 
ciency could be calculated as follows: 

W 2 = ^Re(p 2 • U 2 ), W 5 =±Re(p 5 -U;\ >, = ^ (7) 


4. Simulation results and discussion 


Fig. 3. Flow chart of FPSE with a thermal buffer tube (A. regenerator, B. thermal 
buffer tube, C. displacer, D. power piston; 1. inlet of the regenerator, 2. outlet of the 
regenerator, 3. inlet of the displacer, 4. outlet of the displacer, 5. outlet of the 
engine). 


A modified free piston Stirling engine with hundreds of watts 
acoustic power output is proposed and studied in this article. The 
effect of the size of the thermal buffer tube on the engine 
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performance is investigated. Table 1 shows its dimensions and 
operating parameters. 

We set the length of the thermal buffer tube being 70 mm and 
the diameter of the thermal buffer tube varying from 0 to 50 mm. 
The following figures show the effect of the thermal buffer size on 
the thermodynamic performance. 

Fig. 4 shows the curve of the output acoustic power versus the 
diameter of the thermal buffer tube. As seen in Fig. 4 that as the 
diameter increases, the output acoustic power decreases. The ther¬ 
mal efficiency and relative Carnot efficiency slightly decrease as 
well (Fig. 5). The output phase and volume flow rate increase 
apparently (Figs. 6 and 7). It can be seen that the output acoustic 
power drops slightly although the phase gets worse, which could 
be attributed to the increase of output volume flow rate. However, 
it is adverse to match a linear alternator, so the size of the thermal 
buffer tube should not be too large. 

On the other hand, the displacement of the working gas in the 
thermal buffer tube should not be too large; otherwise, there will 
be a huge heat loss. Define X 2 and X 3 as the amplitude of the work¬ 
ing gas displacement at the inlet and outlet of the thermal buffer 
tube respectively. Their expressions are 


m\ x _ m\ 

ncod 2 ’ 3 ncod 2 



Their values are listed in Table 2. 

Referring to the experiences of a pulse tube refrigerator, the gas 
displacement usually does not exceed 10% of the length of the ther¬ 
mal buffer tube (pulse tube) to avoid large heat loss. Since the 
length of the thermal buffer tube is 70 mm, the displacement 
should be less than 7 mm. As indicated in Table 2, the displace¬ 
ment is too large when the diameter is less than 40 mm. From Figs. 
4-7, it can also be found that the output power and the coupling 
with a linear alternator will be seriously affected under this condi¬ 
tion. Therefore, to avoid large heat loss, the size of the thermal buf¬ 
fer tube should not be too small. 


Table 1 

Structure and operating parameters. 


Structure parameters 


Operating parameters 


Diameter of displacer 

50 mm 

Heating temperature 

923.15 K 

Diameter of rod 

20 mm 

Ambient temperature 

303.15 K 

Diameter of regenerator 

33.5 mm 

Working gas 

helium 

Length of regenerator 

40 mm 

Working pressure 

3 MPa 

Wire diameter 

0.051 mm 

Pressure wave of the 
inlet of regenerator 

0.3 MPa 

Mesh number 

180 

Frequency 

50 Hz 

Porosity 

0.7 

Damping coefficient 

25 kg/s 

Moving mass 

0.185 kg 

Spring stiffness 

15 kN/m 



Diameter of thermal buffer tube /mm 


Fig. 5. Thermal efficiency vs. diameter of the thermal buffer tube. 



Diameter of thermal buffer tube /mm 

Fig. 6. Output phase vs. diameter of the thermal buffer tube. 



Diameter of thermal buffer tube /mm 
Fig. 7. Output volume flow rate vs. diameter of the thermal buffer tube. 



Diameter of thermal buffer tube /mm 


Table 2 

Gas displacement amplitude at inlet and outlet of the thermal buffer tube. 


d (mm) 

X 2 (mm) 

X 3 (mm) 

5 

455.89 

456.82 

10 

113.44 

114.38 

15 

50.04 

50.99 

20 

27.86 

28.82 

25 

17.61 

18.59 

30 

12.06 

13.07 

35 

8.73 

9.76 

40 

6.59 

7.65 

45 

5.15 

6.23 

50 

4.14 

5.24 


Fig. 4. Output acoustic power vs. diameter of the thermal buffer tube. 
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5. Conclusions 

A modified FPSE with a thermal buffer tube is proposed in this 
article. Moreover, a simplified thermoacoustic model is established 
for the new design. Based on the model, the effects of the thermal 
buffet tube size on the performance of the modified FPSE are inves¬ 
tigated in detail. The analysis indicates that the thermal buffer tube 
would not significantly decrease the acoustic power output but 
keep a relatively high efficient operation. Taking into the consider¬ 
ation of several advantages of the modified FPSE, this new config¬ 
uration is worthy of testing and developing in future. 
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